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Abstract 



Motivated by recent experimental studies of Hg and Pb monolayers on Cu(OOl) 
we introduce a zero temperature model of a monolayer adsorbed on a square sub- 
strate. Lennard- Jones potentials are used to describe the interaction between pairs 
of adlayer-adlayer and adlayer-substrate atoms. We study a special case in which 
the monolayer atoms form a perfect square structure and the lattice constant, posi- 
tion and orientation with respect to the substrate can vary to minimize the energy. 
We introduce a rule based on the Farey tree construction to generate systematically 
the most energetically favored phases and use it to calculate the phase diagram in 
this model. 
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1 Introduction 



The study of ordered phases of films of one layer or less on metal substrates is 
an active topic of research in two-dimensional physics [I]. Due to improved film 
deposition and detection techniques, it is now possible to investigate more systems 
and with better resolution than has been possible before f2j. Experimental studies 
of ordered overlayers show a great variety of commensurate and incommensurate 
phases ||, |J. Here we define a commensurate phase as an ordered structure for 
which it is possible to find a primitive cell with basis vectors which are linear 
combinations of lattice vectors of the substrate. If this primitive cell (supercell) is 
larger than the overlayer unit cell and contains several atoms, then we will call this 
a high-order commensurate structure. When the number of atoms in this cell goes 
to infinity, the high-order commensurate structure becomes incommensurate. 

With some exceptions, structures of strongly chemisorbed systems are either 
fully commensurate (every adatom is in registry with the substrate lattice) or in- 
commensurate phases; while high-order commensurate phases are less often seen 
among systems in which strong adatom localization occurs. Because of the lo- 
calized nature of the adsorption and the poor understanding of the interaction 
potentials, theoretical methods to describe the structure and phase boundaries of 
these layers have relied by and large on lattice gas models and the parametrization 
of interaction energies || U . 
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In the opposite limit, when incommensurate structures are formed, other ap- 
proaches have been tried. Calculations by Novaco and McTague and Shiba || 
showed that orientational epitaxy might arise in an incommensurate layer adsorbed 
on a substrate; this effect has been observed for rare-gas monolayers adsorbed on 
graphite and for certain alkali monolayers on metal substrates. 

In this work we introduce a model that describes an ordered layer with square 
symmetry adsorbed on a square substrate. This model is used to compute the zero 
temperature phase diagram for this system. 

2 Experimental Background 

Recently, high-order commensurate structures have been detected for Hg and Pb 
adsorbed on Cu(OOl) |9|, [H]]. In addition to the fully commensurate C{2 x 2) phase, 
two high order commensurate square phases have been observed and studied in 
details: the C(4 x 4) phase (Hg on Cu) and the (5 x 5)i?tan _1 (3/4) phase (Pb 
on Cu). Only one of the atoms of the supercell is in the ideal registry position, 
while all the others sit in lesser energetically favorable locations. The experimental 
results indicate that in these phases the adlayer forms an almost perfect square 
lattice with only small local relaxations caused by the substrate potential. 

Surprisingly, these high-order commensurate phases possess square symmetry. 
We recall that for an unsupported layer with Lennard- Jones pairwise interactions 
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the most stable structure has hexagonal symmetry fll|| . It is then clear that these 
high-order commensurate square structures must be stabilized by substrate forces 
and/or by more complicated interactions, such as three-body forces or angular 
dependent, covalent-like interactions 0. 

3 Classification of Square Phases 

We have developed a system, based on the two dimensional analog of the Farey 



tree of rational numbers [12], to classify all the commensurate square phases on a 
square substrate. Consider a square monolayer on a square substrate such as the 
(5 x 5).Rtan -1 (3/4) phase (Fig. 1). The unit cell of the monolayer is a square, 
which is typically rotated with respect to the substrate main symmetry directions. 
To characterize the structure we pick a side of the square and project it on the 
x and y directions. Both projections are integer multiples of the substrate lattice 
constant, which will be given by the integers p x and p y respectively. The linear size 
of the supercell, in terms of the monolayer lattice constant, is given by q. Therefore, 
one can identify each square phase by the pair of rational numbers 



Px Py 



Using this observation we propose a systematic hierarchical procedure to generate 
a class of commensurate square phases. The advantage of this method is that it 
generates a sequence of phases from the most stable and lowest order ones up to 



the least energetically favorable. 

The list of phases is constructed by starting with a set of four lowest order 
phases (furthest corners in Table 1). Then between each pair of neighboring phases 
one adds a new phase for which the new p x ,Py,Q are the sums of the corresponding 
quantities in the neighboring phases. The four new phases are located at each side 
of the square (midway between corners in Table 1). A new phase is also introduced 
in the center, using as input the two phases on one of the diagonals. Here one should 
choose the diagonal that includes a phase which belongs to a lower generation. 

One can continue this procedure by breaking the new diagram into four separate 
smaller squares and repeating the procedure for each one of them. The main advan- 
tage of this procedure is that it generates phases of low order first (which are most 
likely to be stable) while higher order phases appear in later generations. This con- 
struction should provide a useful guide for the experimentalist willing to explore 
interesting yet undetected structures which are in the vicinity of experimentally 
observed ones. 

4 The Model 

In order to obtain a better understanding of which adlayer structures one should 
expect to find on a square substrate, we have studied a zero temperature model and 
calculated its phase diagram. In particular we consider a two dimensional layer of 
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atoms with Lennard- Jones interactions, subject to a periodic underlying potential 
due to interactions with the substrate. Since the Lennard- Jones interaction is of 
long range, each atom in the overlayer interacts with all the other atoms of the 
overlayer and the substrate. The lattice constant, orientation and position with 
respect to the substrate can all be changed to minimize the energy. Since the 
square phases detected experimentally seem to show little relaxation (and this is 
born out by a mean- field type of calculation ||11||), one can consider that the adlayer 
atoms are constrained to a perfect square structure. 

In this simplified case one can use results for lattice sums to add the adlayer- 
substrate interactions into an effective potential. The energy will then be: 

£ = E V Lj(r m ~ f n ) + E V s(r n ) (2) 

where 



in, ri 



V L j{r) = ^e ad 



is the Lennard- Jones potential between the adlayer atoms, Vs is the underlying 
substrate potential and r n = (x n ,y n ) is the position of the nth adlayer atom. The 
substrate potential is given by: 

V s (f n )=A + J2A l :e^ (4) 

k 

where the sum runs over all the reciprocal lattice vectors of the substrate. The 
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Fourier coefficients are given by 



A = 8n£a 2 

and 



O \ Zq J \ZqJ 



(5) 
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A ' = lOoT {^) K ^ - K ^ (6 > 
where k = \k\, K p / q is a modified Bessel function of order p/q, z is the distance 
between the adlayer and substrate and a s b = (2Ci 2 /C 6 ) 1 / 6 -cr sb is the lattice constant 
of the substrate where C n = ~ l/(^ 2 + j 2 ) n//2 . Here, the parameters a and e for 
the adlayer-substrate interaction are given by cr = (a a d + cr sb )/2 and e = ^e a d e s b- 
For simplicity we choose e s b = 1 and a s b = 1. The distance = a between the 
adlayer and the substrate is then determined so that the zeroth order term, Aq, in 
the potential is minimized. Note that we assume here that the distance z between 
the monolayer and the substrate is the same for all adlayer atoms, although in 
reality these heights vary (0.1-0.2 A) |TJJ and depend on the position of the atom 
with respect to the substrate. 

To study the phase diagram we first choose a structure and then calculate 
the energy per atom. We compare these energies for a set of commensurate and 
incommensurate square phases and obtain the phase diagram for the physically 
interesting regime of 1 < a a d < 2 (Fig. 2). The parameters in the phase diagram 
are a aa i, and the latter being proportional to the strength of the underlying 

substrate potential. We normalize the lattice constant of the adlayer in units of 



the substrate periodicity and define a = a ad /a s b- The phase diagram contains both 
commensurate and incommensurate square phases. Each commensurate phase is 
stable in a tongue of a finite width, in which a = ^jp\ + Py / q, while incommensurate 
phases exist in the gaps between these tongues. In Fig. 2 one can see that the 
tongues start very narrow and as the parameter increases they become 

broader until they intersect in triple points. Above the triple point there is a first 
order transition line between the two tongues, while all the incommensurate phases 
between them disappear. The phase diagram is asymmetric and this is due to the 
hard core nature of the Lennard- Jones potential. The adlayer phase can expand 
considerably in order to gain energy from the underlying substrate potential but 
it cannot contract as much due to the hard core repulsion. Also note that the 
high order tongues are extremely thin. This is due to the asymptotic properties of 
the modified Bessel functions, which appear in the Fourier coefficients A\.. For high 
order harmonics of the potential the amplitudes decay exponentially, or oc e~ z ° k . 
Thus, the energy gain of high order commensurate phases becomes rapidly very 
small as the order of commensuration increases. 

Finally, we have calculated the phase diagram with the inclusion of (incommen- 
surate) hexagonal phases; in this case, most of the square phases are washed out, 
except for the fully commensurate ones. This is to be expected, since high order 
phases appear for weak interactions with the substrate (lower part of Fig. 2) and 



8 



the energy gain for commensuration is small. 

5 Conclusion 

In summary, we have studied a simple model for the energetics of square monolayer 
phases on a square substrate. To study the phase diagram we introduced a simple 
classification that allowed us to find the most energetically favored phases. We find 
that when only square phases are considered, each one of them has a tongue in which 
it is stable. The C(4 x 4) phase (Hg on Cu) and the (5 x 5)i?tan _1 (3/4) phase 
(Pb on Cu) belong to very narrow tongues (not shown in Fig.2) which intersect 
the a a d axis at 1.265 and 1.25, respectively. The fact that hexagonal phases win 
out (in our calculation) might be due to a variety of factors, including the fact 
that the observed phases are metastable and the actual interactions might be more 
complicated than the ones used in this study. 

One of us (G.V.) would like to acknowledge partial support from NSF grant 
8802512. 
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Figure Captions 

Figure 1: The (5 x 5)i?tan _1 (3/4) phase that appears in Pb on Cu(OOl). Empty 
circles represent substrate atoms, while full circles represent adlayer atoms. 

Figure 2: Phase diagram (at K) of two dimensional homogeneous square phases 
on a square substrate, with Lennard- Jones interactions. Note the clear asymmet- 
ric shape, which is due to the asymmetry of the Lennard- Jones potential. It in- 
cludes three fully commensurate square phases: 1 x 1 (or (0/1,1/1)), C(2 x 2) 
(or (1/1, 1/1)) and 2 x 2 (or (0/1,2/1)). It also includes a few examples of the in- 
finitely many high order commensurate square phases (from left to right): [1/2, 2/2], 
{3/5,6/5}, [0/2,3/2], [1/2,3/2], and {3/4,7/4}. 

Table Caption 

Table 1: A classification for square monolayer phases on a square substrate. In this 
hierarchical classification, the lowest order phases, which are typically the most sta- 
ble, appear first, while higher order phases appear in later generations. Generations 
are distinguished by different symbols: (),[],(), and { } for generations 1, 2, 3 
and 4, respectively. For example, the experimentally observed high order commen- 
surate square phases C(4 x 4) and (5 x 5)i?tan^ 1 (3/4) appear here as {2/5,6/5} 
and {3/4,4/4} respectively. 
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